Quantum Kinetic Theory of Current-Induced Torques in Rashba Ferromagnets 
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Motivated by recent experimental studies of thin-film devices containing a single ferromagnetic 
layer, we develop a quantum kinetic theory of current-induced magnetic torques in Rashba-model 
ferromagnets. We find that current-induced spin-densities that are responsible for the switching 
behavior are due most essentially to spin-dependent quasiparticle lifetimes and derive analytic ex- 
pressions for relevant limits of a simple model. Quantitative model parameter estimates suggest 
that spin-orbit coupling in the adjacent metal normal magnetic layer plays an essential role in the 
strength of the switching effect. 
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Introduction — The central goal of spintronics is the dis- 
covery of efficient mechanisms for electrical control of 
nanomagnet orientation. The past fifteen years have 
witnessed many advances (see reviews [IH3], and refer- 
ences therein) that are potentially important for infor- 
mation storage technologies and, because they depend 
on improved understanding of non-equilibrium collective 
properties, also interesting from a fundamental point of 
view. A breakthrough has recently been achieved with 
the demonstration [H [5] of reliable electrically-controlled 
magnetization switching in trilayer thin- film devices, il- 
lustrated schematically in Fig. [I] that have a single 
perpendicular-anisotropy magnetic layer. (Magnetiza- 
tion control with weak current induced torques had been 
observed previously[6 in (Ga,Mn)As.) In the demonstra- 
tion devices the ferromagnetic layer was Co, the normal 
metal layer Pt, and the insulator AI2O3. (Ref. [5] also 
studied other insulators to partially demonstrate the uni- 
versality of the effect.) 

Although they appear to be observing the same effect, 
Refs. [4] and [5] ascribe the current induced torques to 
different mechanisms. Ref. [4 speculates that the effect 
is due to Rashba spin-orbit interactions of carriers in the 
ferromagnet and that it is analogous to the surprising 
semiconductor two-dimensional electron gases (2DEGs) 
z direction spin-densities induced by current flow parallel 
to an applied magnetic field [TJ [8;. On the other hand, 
Ref. [5 argued that the torques are due to a large spin 
Hall current in Pt that flows vertically into the magnetic 
layer. The two effects are identical from the symmetry 
point of view, as already observed in Ref. [4 . 

The spin-Hall-effect point of view has been elaborated 
on rather extensively in Ref. [5]. The model we study in 
this paper is therefore most closely related to the Rashba 
effect interpretation. In general, current-induced torques 
in ferromagnets have for several years been recognized 
theoretically as a general consequence of spin-orbit cou- 
pling (see Refs. [9HTT], as well as review [3 , and references 
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FIG. 1: (color online) Schematic view of trilayer devices for 
electrical magnetic switching. A perpendicular-anisotropy 
ferromagnetic layer (Co) is sandwiched between a metallic 
layer (for example Pt) and an insulating ( for example AI2O3) 
capping layer. The transport current, jt r , flows in the itiner- 
ant ferromagnet and normal metal in the direction parallel or 
antiparallel to the external magnetic field, B e *t- 



therein). For Rashba ferromagnets with z direction inver- 
sion symmetry breaking, out-of plane spin-densities can 
be explained [8 only by a full solution of the two-band 
quantum kinetic equation. With this motivation we re- 
port on a full quantum kinetic theory of current-induced 
torques in a ferromagnetic conductor with structural in- 
version asymmetry. 

Rashba Ferromagnets — Quantitative microscopic theo- 
ries of transport in transition metal ferromagnets are 
complicated by the presence at the Fermi energy of sev- 
eral 3d- and 4s-derived electronic bands that have com- 
plicated Fermi surfaces. Most transport properties of 
collinear ferromagnetic conductors, including the impor- 
tant giant and tunnel magnetoresistance effects, can nev- 
ertheless be understood qualitatively in terms of a simple 
two-channel model [T3UT6] in which ^-spin and ^-spin car- 
riers conduct in parallel with conductances that depend 
strongly on the magnetization direction and its space de- 
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pendence. The direction dependence is due to a combi- 
nation of spin-split band structure and spin-dependent 
scattering effects. We will see that the spin-dependent 
transport is also important for the present problem. We 
take the point of view that qualitative insights can be 
achieved by studying a simplified model that accounts 
for the interplay between structural inversion asymme- 
try and spin-dependent transport. 

Cobalt is a strong ferromagnet [17 in which the 
majority-spin Fermi surface intersects only a high- 
velocity free-electron-like (s) band with a substantial 
conductance, while the minority-spin Fermi surface in- 
tersects several complex bands with generally small ve- 
locities, large densities-of-states, and dominant atomic- 
like d-character. Minority- spin states at the Fermi en- 
ergy have strongly hybridized free-electron-like and d am- 
plitudes. Even high- velocity minority- spin states con- 
tribute less effectively to transport because they scat- 
ter rapidly into the many available d-states. To cap- 
ture the essential physics we employ a phenomenolog- 
ical two-band model with short-lifetime minority spins 
and long-lifetime majority spins, neglecting the d-orbital 
contribution to transport as is customary in theories of 
GMR [T3HT6] . Scattering is spin-dependent because any 
crystal defect will produce different potentials for major- 
ity and minority electrons when strong exchange is taken 
into account, and because minority s-electrons can scat- 
ter into empty minority d-electron states. As is common 
in model studies, we add spin-orbit coupling and vertical 
structural inversion asymmetry to the two-band model 
by adding a Rashba SO term to the Hamiltonian. 

These considerations lead to a Rashba-model ferro- 
magnet with band Hamiltonian: 
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Here e p denotes the operator for the band energy of spin- 
orbit and exchange uncoupled itinerant electrons. Be- 
cause the Co layer is thick enough (~ 10A) to support 
several 2D subbands, we take p to be a 3D vector. The 
momentum-dependent effective magnetic field includes 
an exchange contribution in the direction m of the d- 
electron magnetization and a spin-orbit field in the z x p 
direction: B p = — A xc m + 2ai x p, a being the strength 
of the Rashba coupling. (We neglect Zeeman coupling in 
comparison, assuming that it is important only in help- 
ing to fix the direction of m.) To avoid confusion in what 
follows, we define m to be the direction of local magneti- 
zation. Thus the majority (minority) electrons have their 
spins antialigned (aligned) with the local magnetization, 
which corresponds to A xc > 0. Udi s is the disorder po- 
tential. We adopt a short-range-disorder model with a 
spin-dependence that is dictated by m: 



(2) 



where the index i labels the impurity positions r^, P± = 
(1 ±6- m)/2 are projectors onto the local direction of the 
exchange field, and and u± characterize the strengths 
of majority and minority spin scattering, respectively. 
Because majority-spin electrons have a higher mobility 
we must chose u± > u^. Our model neglects inelastic 
phonon and magnon scattering which should be relatively 
unimportant in highly disordered ultrathin films even at 
room temperature. Below we refer to the high- velocity 
electrons described by the Rashba ferromagnet model as 
the transport electrons, and to the spin-density of the 
d-electrons as the magnetization. 

Current Induced Torques — Because the dynamics of the 
itinerant electrons in magnetic conductors is much faster 
than that of the magnetization, it is sufficient to consider 
the kinetics of transport electrons in a static exchange 
field. 

The effective magnetic field exerted on the magnetiza- 
tion by the transport electrons is due to their mutual ex- 
change interaction [10] and therefore proportional to A xc : 



eff 



1 S(H) S 
M s Sm 



M s 



(S)s 



(3) 



The subscript V on the angle brackets indicates the 
expected value in the transport steady state. Here 
s is the transport electron spin-density operator s = 
Mj.B crcr /^ r cr //2 and Ms is the d-electron spin-density. 
The theory of current induced torques is therefore equiv- 
alent to a theory of (s) s . The required transport theory 
differs essentially from the 2DEG problem considered in 
Ref . [8] , since the spin dependence of transport comes pri- 
marily from the spin-dependence of impurity scattering 
in an itinerant ferromagnet and the presence of a minor- 
ity Fermi surface for d-like electronic bands, rather than 
solely from band-structure spin-splitting and angular de- 
pendence of scattering as in a 2DEG [8]. (Below we will 
also present the leading-order results for the latter case.) 
The component of (s) s perpendicular to the magnetiza- 
tion, which is responsible for magnetization switching in 
the devices of interest, depends essentially on current- 
induced interband coherence and must be evaluated us- 
ing a two-band quantum kinetic theory. 

The kinetic equation can be obtained along the route 
outlined in Ref. [12]. In this paper we do not, however, 
switch to the energy representation, choosing to work in 
the more intuitive p-representation. The final form of 
the equations we use is therefore slightly different. 

To facilitate the discussion that follows, we first es- 
tablish some notation. In the absence of disorder, the 
eigenvalues of Hamiltonian ([!]) are given by 
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where the index v = ± distinquishes majority (+) and 
minority (— ) bands. The projection operators onto the 
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corresponding eigenstates are 



where the check accent denotes matrices in Keldysh 
space [18] . 



1 ( 6-Bp 



(5) 



For weak disorder of the form specified in Eq. ([I]), the 
disorder self-energy is 

S(r) = n,7G(r,r)7, 7 = + u t P°, (6) 



The quantum kinetic equation for the two-band density 
matrix / p , linearized with respect to a time-independent 
uniform external electric field E, is 



dtf P + \ - ^pBpfi, d r f p J - ^[B p B, / p ] + eEa p /; q = I st . (7) 

Here / p is 2 x 2 matrix in spin-space and — S ^ JV P pv fth(tpv) , where fth(E) is the Fermi function, is the equilibrium 
distribution function. Disorder effects appear in the collision integral I st which takes the form 

I at = -TTTli j ^~)3 S ^ 6 P V/ ) (PpvfpjPp'v'J + iPp'v'lfpPpv ~ Ppvlfp'Pp'vH jPp'v* /p'T^pi/) • 



It is easy to check that the collision integral vanishes 
when f p is replaced by /p q . 

In what follows, we assume that the spin splitting of 
the transport electrons B p is small compared to their 
Fermi energy Ep. In the analogous semiconductor 2DEG 
kinetic equation [8] one has to keep at least the terms lin- 
ear in B p /Ep in the kinetic equation in order to con- 
sistently describe the coupling between charge currents 
and spin densities. In the present case, however, the spin 
and charge densities are already coupled at (B p /Ef)° 
level because of spin-dependent scattering. In order to 
isolate the essential physics as transparently as possible 
we initially neglect Fermi surface splitting in the field- 
generation and collision terms. It is legal to keep the 
precession term, which is linear in B p , since it competes 
with relaxation rates only. We comment further on the 
influence of spin-split bands below. 

We now specialize to the stationary and uniform case. 
After straightforward manipulations, using P pu = 1, 
performing the integrals over p' in the collision integral, 
and representing the distribution function as a sum of 
scalar and vector parts using f p = n p + fi • f p , we obtain 
the following equations: 

1 _ HI 

eEd p fth = (n p -n p ) (f p - f p ), 

m 1 

5 p xfp = (n p - n p ) (f p - f p ) 

m x f p x m, (8) 



where l/r s ^ = 7rriiNo{u^ ± ui) > 0, Nq being the den- 
sity of states at the Fermi level, and the overbar accent 
denotes the average over the directions of p. The spin- 
decoherence rate, 1/tj_, i. e. the rate of local relax- 
ation of spin-components perpendicular to m, is equal 
to 7rriiNo(u_i — u^) 2 in this model, and thus is a func- 
tion of l/r Sj d. We keep it as an independent parameter 
to recognize that this property does not hold for general 
spin-dependent disorder models [22]. Longitudinal spin 
density (aligned with m) is not generated to linear order 
in the spin-orbit interaction in the present case and we 
therefore do not need to the introduce the corresponding 
relaxation time. 

The solution of Eq. ([8| is straightforward. The scalar 
equation is used to eliminate the charge response from the 
vector equation, thereby introducing a generation term 
and an additional relaxation term in the equation for 
f p . To linear order in a, we first obtain an equation for 
the longitudinal spin response (the component of f p par- 
allel to m), which does not directly produce a torque. 
When this solution is substituted into the equation for 
the transverse response it adds a generation term to the 
equation for the perpendicular component because of its 
precession in the spin-orbit field. The current-induced 
spin-densities and currents are then obtained by sum- 
ming the scalar and vector responses over p. 

Using Eq. (|3|, we finally obtain the following expres- 
sions for the effective magnetic fields that act on the fer- 
romagnet's magnetization, and for the current density 
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flowing the in the film: 
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where j tT is the transport current density, n tr is the trans- 
port electron density per spin; we assumed for simplicity 
a parabolic spectrum of transport electrons, and intro- 
duced m, the correpsonding band mass (the same for 
majority and minority electrons). The latter assumption 
is not crucial, and has been invoked only for reasons of 
clarity of the expressions. Finally, rVl = r s _1 ± r^ 1 are 
the minority (1/rjJ and majority (1/r-f) spin scattering 
rates. These expression are the main result of this work. 
We identify the field H^, which produces a torque iden- 
tical to a field in the y direction and cannot switch a 
perpendicular film as the Rashba field. We note that 
the expression for the Rashba field given here is para- 
metrically larger than the one in existing literature [9], 
which contains an additional A KC /Ep smallness factor. 
If Ref. [19] was analyzed using our expression it would 
decrease the experimental estimate of the Rashba cou- 
pling parameter. The field Hg ff has the same dependence 
on magnetization as the effective field produced by a y- 
polarized spin-current flowing into the magnetic layer, 
and can switch perpendicular currents [4j [5] . 

In the general case of spin-dependent scattering and 
exchange fields of arbitrary strength the two-band quan- 
tum kinetic equation cannot be solved analytically. 
There are however other limits in which instructive an- 
alytic results can be obtained. One interesting limit is 
that of a 2D Rashba ferromagnet with spin-independent 
disorder and a spin-splitting that is larger than the Bloch 
state lifetime. Results for the effective field, h e ff, to the 
leading order in a and in the clean limit, B p r ^> 1, are 
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In the above expressions 70 = (^f^; — is the 
nonparabolicity parameter, introduced in Ref. |8]. In 



Eqs. (10), the field is the Rashba field found earlier 



are proportional to 70 but even in m z since they can- 
not lead to switching. Terms with the symmetry of H^ ff 
from Eq. ^ do appear at 0(a 3 ) order, and are also pro- 
portional to 70, but the corresponding expressions are 
cumbersome and are not shown. 

By examining these two analytically accessible lim- 
its, we have concluded that the crucial ingredient in the 
switching behavior is spin-dependent scattering. It seems 
likely to us that the most important source of this spin- 
dependence in Co is the presence at the Fermi energy of 
minority spin d electrons which is captured in Eqs. [9] 
Discussion - The physical origin of the effects described 
here is the fact that when an electric field is applied to 
a ferromagnet with different majority and minority mo- 
bilities, the two spin species react to it differently. This 
creates spin polarization in momentum space, whose an- 
gular dependence "matches" the angular dependence of 
the y component of Rashba SO field. Thus a compo- 
nent of spin polarization perpendicular to the ferromag- 
net's magnetization is generated. This component then 
precesses around the exchange field, and gets relaxed ei- 
ther by magnetic scattering or the Dyakonov-Perel mech- 
anism. A similar scenario was outlined previously [8] in 
connection with the explanation for out-of-plane spin po- 
larizations in in-plane fields in 2DEGs with Rashba in- 
teractions. 

Despite the similar underlying mechanism, the leading 
order results in the two models we considered, Eqs. ([9| 
and (10) have different form. In particular, Eqs. (l9|) cor- 



respond to the case of spin-dependent impurity scatter- 
ing, and Eqs. (10) are those for the case of a spin-split 



band structure and spin-independent impurity scatter- 
ing. The difference is thus not surprising: since the ori- 
gin of the spin-dependent scattering is different in the 
two models, generation and relaxation terms have differ- 
ent dependence on the Rashba SO strength. 

There is some correspondence between the current- 
induced fields and torques that we have evaluated for 
a Rashba ferromagnet model, and those found experi- 
mentally. The torques exerted on the magnetization by 



H eff , Eqs. (|9j), have the same symmetry properties as the 
experimental current-induced torques in Refs.[4] and [5]. 
In addition, the field has the same sign as observed ex- 
perimentally, since it implies switching from up to down 
direction for current flowing in the direction of the ap- 
plied external in-plain magnetic field, if a is positive. 
The fact that it is positive in the experiments of Miron 
et al. follows from their observation that points in 
the positive y direction (the coordinate system used in 
this work is identical to theirs). This universality is en- 
sured by the transport properties of Co, in which the 
majority electrons have a much larger lifetime. In the 



in Ref. [9 . The field h^ ff , despite being aligned with the 
magnetization, can in principle lead to switching since 
it is odd in m z . The factor 1 + m 2 z in the denominator 
of the expression for h^ ff can be traced to the fact that case of a 2DEG, Eqs. (10), the direction of the switching 
for a 2DEG with a Rashba SO, the Dyakonov-Perel [20] 
in-plane spin-relaxation time is twice as long as the out- 
of-plane one. We have dropped contributions to h e ff that 



field is less universal, in the sense that it depends on the 
details of the band structure through the nonparabolicity 
parameter 70. 
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The magnitudes of the effective fields in a Rashba fer- 
romagnet, Eqs. (|9|, depend on two parameters, a and 
Axc^x, which could be estimated by detailed first prin- 
ciples studies of Co/Pt multilayers in the presence of a 
structural asymmetry. By comparing the ratio of switch- 
ing and Rashba fields in Eqs. ([9|, we observe that the 
values of Hf ff /j tr w lO^TWA- 1 claimed in Ref. [19], 
and Hf ff /.7tr ~ 10- 9 Tcm 2 A-\ claimed in Ref. g] in a 
similar setup, yield a plausible value of A xc t_l ~ 10. 

For a large A xc t_l our results for the Rashba spin- 
density correspond to 

s R ap F v D 

^ (11J 

n tr A xc v F 

where vd ~ eEr/m is the transport drift velocity. In 
Co we estimate that vjj ~ 100 m/s at a current density 
of 10 12 A/m 2 , compared to a transport electron Fermi 
velocity ~ 10 6 m/s. To estimate that Rashba coupling 
parameter of Co on Pt we note that the majority-spin 
Fermi level is ~ O.beV above the top of the d-bands in 
Co and and below the top of the band in Pt. This im- 
plies an electric potential drop comparable to the Fermi 
energy over the thickness tco of the Co layer and a 
Rashba coupling strength at the Fermi energy ap F that 
is ~ E F /(mc 2 p F t Co) which should be around 10 5 times 
smaller than A xc . We conclude that it is hard to ex- 
plain dimensionless current-induced spin-densities much 
larger than ~ 10 -9 at 10 12 A/m 2 in Co based on this bulk 
mechanism. On the other hand the quantitative analysis 
in Ref. [5 suggests a switching dimensionless spin-density 
at this current level that is ~ 10 -6 . We thus note that 
a substantial enhancement of SO interactions is needed 
to explain the experimental switching current densities 
based on the Co Rashba mechanism. This enhancement 
could come from the hybridization of Co transport elec- 
trons with Pt d-orbitals, analogous to SO enhancement 
by surface alloying with heavy elements [21]. This hy- 
bridization must play an essential role in the switching 
effect in both Rashba coupling and spin-Hall scenarios. 
We therefore propose that the effect is strongest when 
these bands are closely aligned. If so, this may explain 
the efficacy of Co/Pt structures and also suggest strate- 
gies for optimizing the effect. 
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